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Cosmological Simulation of 
Gravitational Weak Lensing in Open Models 

Kenji Tomita 

'^q' ■ Yukawa Institute for Theoretical Physics, Kyoto University, Kyoto 606-8502 

^\ , Cosmological simulation of gravitational weak lensing in open models with (J7o,Ao) = 

(0.2, 0) and (0.4, 0) is treated using the direct integration method, which is independent of the 
multi-lens-plane method. Statistical averages of optical quantities (the convergence, shear, 
and amplification) are derived for ray bundles passing through the clumpy universe, in which 
the initial state has the CDM spectrum with the power n = 1, and all particles are regarded 
as equivalent (compact) lens objects with a galactic size. The compact lens objects consist 
of galaxies and dark matter particles. The evolution of the spatial mass distribution is given 
by A^-body simulations in a tree-code. This paper is an extension of our previous one in flat 
models. The behaviors of optical quantities are shown in comparison with a low-density flat 
model with (0.2,0.8). It is found as a result that (1) optical quantities in the open model 
increase with Qq, and (2) they are larger than those in the flat model with the same i?o, 
and the difference is comparable with that between two open models with Qo = 0.2 and 0.4. 
The comparison of the optical quantities with the corresponding observed ones will give us 
some information on the observational value of f2o. 
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^ ■ SI. Introduction 
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Gravitational weak lensing is one of the most useful tools for probing the inho- 



Q^. mogeneous distribution of mass in the universe. The small deformation in the shape 

Q I of distant galaxies caused by gravitational tidal force is directly connected with the 

^ ' total mass distribution in the neighbourhood of light paths to galaxies. Recently_thB- 

c/3 . oretical and observational studies on the deformed images of distant galaxies "13 

^ ' have shown that in a near future the lensing information on the mass distribution 

brought with large telescopes will be useful to put severe constraints on dark matter, 
^ . dark halos and allowed ^alues of cosmological parameters. 

^1 In a previous papero* we treated the statistics of gravitational weak lensing in flat 

cosmological models in the direct integration method, in which the ray shooting was 
performed by solving null-geodesic equations numerically in the intervals between an 
observer and sources. Here weak lensing does not mean the weakest limit of lensing, 
but means that we treat only rays without caustics. 

In this paper we extend the treatment to the case of open models and derive the 
average values of optical quantities on various angular scales, in comparison with 
those in the flat case. The largest difference of treatments in curved cases from 
the flat cases is that in curved cases we can construct no model universes covered 
everywhere with periodic boxes, and so, for the calculation of the gravitational forces 
from particles, we must assume a set of periodic boxes which are closely connected 
only along each light ray. 

In §2 we show treated model universes and the equations for ray shooting in 
open models, and in §3 derive the statistical averages of optical quantities. Their 
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behaviors are shown in connection with those of the angular diameter distances. In 
§4 the comparison of average optical quantities with observed ones is discussed for 
the estimate of i7o- In Appendices A the derivation of null-geodesic equations is 
shown in open models. 

§2. Model universes and the ray shooting 

The background model universes are assumed to be spatially open or have the 
negative curvature, and the behavior of light rays passing through inhomogeneities 
in them is considered at the stage from an initial time ti (with redshift zi = 5) to 
the present time to. In the Newtonian approximation, the line-element is expressed 
as 



-(1 + 2^/c^)c^dt^ + (1 - 2^/c^)a'^{t){dxfl[l + K-{xff, (2-1) 



where K is the signature of spatial curvature (±1, 0). In the following we take 
K = —1. The normalized scale factor S = a{t)/a{to) satisfies 



f2Q - {Qq + Ao - 1)S + Ao^^ 



where r = Hgt and ao(= a{to)) is specified by a relation {cHq ^ /a^f = 1 
The gravitational potential if is described by the Poisson equation 



(2-2) 
Qq — Aq. 



'-A^ = [1 



X 



,2i2 



+ 



dx^ ' [1 - l{x)^]^ dx\ 



= A'eGpb\p{x)/pb - 1], 
where Pb{= Pbo/S^) is the background density and 



SttG 

lOO/iMpc"^ km s"^. 



2.77 X 10^^ f2oh^MQ Mpc 



-3 



(2-3) 



(2-4) 



where Hq = 100/iMpc~ km s~ . In our treatment the inhomogeneities arc locally 
periodic in the sense that the physical situation at x is the same as that at a; -|- In, 
where the components of n(= (n^,n^,n^)) are integers. In an arbitrary periodic 
box with coordinate volume l^, there are N particles with the same mass m. It is 
assumed that the force at an arbitrary point is the sum of forces from N particles in 
the box whose center is the point in question, and that the forces from outside the 
box can be neglected. 

In this paper we take two open models (01 model and 02 model) with {f2o, Aq) = 
(0.2,0) and (0.4,0), respectively, and a flat model (L model) with (0.2,0.8) for com- 
parison. The present lengths of the boxes arc 

Lo = a{to)l = 50h-^ , 39.7h-^Mpc (2-5) 

for f2o = 0.2, 0.4, respectively. The particle number is 32^ in both models, and so 

m{= pbqLq^/N) = 2.11 X lO^^/i-^M©. (2-6) 
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The distributions of particles in these jpodels were derived by the numerical 
A^-body simulations using Suto's tree-codeB^ during the time interval z = and 
z = zi, where zi = 5.0 for all models. All particles in these models are regarded as 
equivalent (compact) lens objects with a galactic size, which consist of galaxies and 
dark matter. 

For the compact lens objects we assume the physical softening radius a{t)xs = 
20/i-^kpc. 

Light propagation is described by solving the null geodesic equation with the null 
condition. Here let us use T = ^ ln[a{t) / a{ti)] as a time variable and Tq = {T)t=to- 
Then we have dS = 2exp[2(r - To)]dT, so that 



where 
and 



cdt = RcR [Qq + (1 - - Ao)5 + XoS^]-^/'^e^'^dT, 

R^Lo/[{l + zi)N'/^] 
CR^2{c/Ho)/[R{l + zif/% 



The line-element is 



ds^/R^ = -CR^[Qo + (1 - I?o - Xo)S + XoS^]-^e^^{l + a^)dT^ 
+ (1 - act>)e^^dy^/F{y)\ 

where y° = T, = a{ti)x'/R, ip = {Gm/R)(l), Rq = Ra^/ai = (1 + zi)R, 



2Gm 3 a 



a 







C^R TT (C^)2 



and 



F = l- -{RoHo/cf{l -no- Xo){yf. 



The equations for light rays to be solved are 



dT 



(2-7) 

(2-8) 
(2-9) 

(2-10) 

(2-11) 
(2-12) 

(2-13) 



^ = -[3Aoe^(^-^») + (1 - - Ao)]e2(^-^o)^YG(r) + a^K^ 
dT oT 



-1 T 

ICR e 



d<P/dy'/G{T) - 2g^K^K^ 



+ {RoHo/cfil -f2o- Xo)F'^CRe^ 



+ -(l + 2a(A)/G(r) 



where 7 = a(cR)^ and 



G{T) = Qo + {l-no- Ao)e2(^-^o) + Aoe6(^-^°). 



(2-14) 



(2-15) 
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The null condition is 

^(A'*)2 = 1 + 2q(/>. (2-16) 

i 

The derivation of these equations is given in Appendix A. 

The potential (j) is given as a solution of the Poisson equation. Because the 
ratio of the second term to the first term in the right-hand side of Eq. (2-3) is 
{RQHQ/cY{yY[Ay /\y\\ « 1 for z < zi, the Poisson equation in a box can be 
approximately expressed as 

Fiy^f = 4^GpB[p(^)/PB - 1], (2-17) 

where we used that F{y) ~ 1 for y ~ the box size and can be approximately replaced 
by the central value F{y^ for y » the box size and the suffix c denotes the central 
value in the box. For point sources with p = mJ2n ^{(^^iu ~ Hn)) is the particle 
number), we have (p = 'Pi + 4*2, where 



nyc)e-'^T.u7^^ (2-18) 



\y-y 

and 02 displays the contribution from the homogeneous background density. Here 
let us use for y another coordinates y expressing the space in a locally flat way, 
where the two coordinates are connected as 

y= dy/F{y), (2-19) 
JO 

and the lengths between two points in boxes in two coordinates are approximately 
related as 

Ay = Ay/F{y,). (2-20) 
Then (pi is expressed in terms of y in the usual manner as 

</>, = -e-2^^-J— . (2-21) 

n \y-yn\ 

Corresponding forces are expressed as fi = dcpi/dy^ = [d(pi/dy^]/F{y^). It should 
be noted that the contribution of d(p/dT is negligibly small, compared with that of 

In flat models the universe is everywhere covered with periodic boxes continu- 
ously connected as in Fig. 1. In open models it cannot be covered in a similar way, 
but we can consider only a set of local periodic boxes connected along each light 
ray, as in Fig. 2. In these boxes we can describe the evolution in the distribution 
of particles in terms of local flat coordinates y, because the size of boxes is much 
smaller than the curvature radius. 

The time evolution in the distribution of particles was derived by performing 
the A^-body simulation in the tree-code provided by Suto. The initial particle distri- 
butions were derived using Bertschinger's software COSMICSu under the condition 
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Fig. 1. Light rays and periodic boxes in the flat space. 




Fig. 2. Light rays and local periodic boxes shown schematically in a curved space. 

that their perturbations are given as random fields with the spectrum of cold dark 
matter, their power n is 1, and their normalization is specified as the dispersion 
as = 0.94 with the Hubble constant h = 0.7. 

For the integration of the above null-geodesic equations, we calculate the po- 
tential at a finite number of points on the ray which are given at each time step 
(Z\r). Then particles near one of the points on the ray have a stronger influence 
upon the potential than particles far from any points. To avoid this unbalance in 
the calculation of the potential, we take an average of the potential cpi by integrating 
it analytically over the interval between one of the points and the next point. The 
expression for an averaged potential (p was given in the previous paper. Moreover, to 
take into account the finite particle size as galaxies or clouds, we modify the above 
potential for point sources using the softening radii = a{ti)xs/R- The modified 
potential is produced by replacing (y — to (y — y^)^ -|- {Vsf' in the potential for 
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point sources. 

The initial values of are given so as to satisfy Eq. ( 2-16 ). The integration of 
Eqs. ( |2T3D and (2-14) with the modified potential was performed using the Adams 
method as in our previous papers. As the time step we assumed AT = [ln6/2]/A'^5 
with /Ns = 3000 (in most cases) -10000. 

§3. Statistical behavior of optical quantities 

We treat the deformation of ray bundles over the interval from z = to z = 5 
measured by an observer in a periodic box. In the same way as in the previous 
paper q), we consider here the ray bundles reaching the observer in a regular form 
such that the rays are put in the same separation angle 9, and calculate the change 
in the angular positions of the rays increasing with redshift in the past direction. 
From this change we find the behavior of optical quantities. 

Basic ray bundles consist of 5 x 5 rays which are put in the square form with the 
same separation angle 9 = 2 — 360 arcsec. Many bundles coming from all directions 
in the sky are considered. Here we take 200 bundles coming from randomly chosen 
directions for each separation angle. In order to express angular positions of rays, 
we use two orthogonal vectors e^^^ and e^g) plane perpendicular to the first 

background ray vector {K^)b- Then the angular coordinates [X(m,n),Y{m,n)] of 
25 rays relative to the first ray with (m,n) = (1, 1) at any epoch are defined by 

X{m,n) =J2[y\m,n) -y\l,l)] eJi)/yB(l, 1) + X(l, 1), 

i 

Yim,n)=J2[y\m,n)-y\l,l)] ej2)/yB(l, 1) + y(l, 1), (3-1) 

i 

where = + and ys = [J^iivh)'^]- Since all angular intervals of rays at the 
observer's point are the same (= 9), differentiation of angular coordinates of the 
rays at any epoch with respect to those at observer's points is given by the following 
differences: 

Aii{in, n) = [X{m + 1, n) - X{m, n)\/9, 
^i2(m, n) = [x{m, n + 1) - X{m, n)]/9, 
A2i{m^ n) = [y(m + 1, n) — Y{m, n)]/9, 
^22("^, n) = [Y{m, n + 1) - Y{m, n)]/9, (3-2) 

where m and n run from 1 to 5. From the matrix Aij{m,n) we derive the opti- 
cal quantities in the standard manner, B as the convergence (K(m,n)), the shear 
(7j(m,n), i = 1,2), and the amplification {ii{m,n)) defined by 

K(m, n) = 1 — tr(m, n) /2, 7i(m., n) = [A22{m, n) — ^ii(m, n)]/2, 
72(m,n) = -[Ai2(m,n) + ^2i("i, n)]/2, 

7^ = (71)^ + (72)^ = [tr(m,n)]2 -det(^ij-(m,n)), 
/x(m, n) = 1/ det(Ajj(?Ti, n)), (3-3) 



The average optical quantities in each bundle are defined as the averages of optical 
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quantities for all rays in the bundle as follows: 



/42, (3-4) 



and so on. In this averaging process the contributions from smaller scales can be 
cancelled and smoothed-out. The above optical quantities at the separation angle 6 
are accordingly derived in the coarse-graining on this smoothing scale. 

For the present statistical analysis we excluded the caustic cases and considered 
only the cases of weak lensing in the sense of no caustics. The averaging for all 
non-caustic ray bundles is denoted using <> as < >, < 7^ > and < {fi — 1)^ >. 
Because K{m,n), 'ji{m,n) and fi{m,n) — 1 take positive and negative values with 
almost equal frequency, < k >, < "fi > and < /x > — 1 are small. 

In Figs. 3 and 4, we show the behavior of < > and < 7^ > for various 
separation angles 6 = 2 arcsec — 360 arcsec (= 6 arcmin). It is found that < >i/2 
and < 7^^ >i/2 for [2q = 0.4 (model 02) are larger than those for i7o = 0.2 (model 
01) at all separation angles = 2 — 360 arcsec and the difference A is about 0.005 — 
0.01 at both periods z = 1 and 2, and that for f2Q = 0.2 those in the open model (01) 
are larger than those in the flat model (L) and their difference A' is comparable with 
the above A. The amount of optical quantities for 9 ~ 350 arcsec is consistent with 
the results of Bernardeau et al. (cf. their Figs. 3 and 4)B' and NakamuraB^ in their 
treatments. The differences of quantities between various cosmological models are 
closely connected with the differences of angular diameter distances in their models, 
because for the same separation angle larger angular diameter distances correspond 
to longer inhomogeneities with smaller amounts and so smaller deformations are 
brought to rays passing through them. Here let us examine the behavior of angular 
diameter distances to explain this situation. In weak lensing the Friedmann angular 
diameter distances can be used approximately to most rays, and their behaviors are 
shown in Fig.5 for models with {f2o,Xo) = (1, 0), (0.4, 0), (0.2, 0) and (0.2,0.8). The 
basic equations for deriving this distances can be seen in a standard text.Q^ From 
this figure we find that these models stand in inverse order of length of the distances, 
and so that in model L the distance is larger for equal redshifts and rays with the 
same separation angles have smaller deformations than in model 01. 

In Tables I, II and III, we show for models 01, 02 and L the numerical values 
of < K >, < >^/^, < /U >, < (/i - 1)^ >^^^, < 71 > and < 7^ >^/^ at epochs 
z = 1,2, .., 5 for 6* = 2 arcsec. 



§4. Concluding remarks 



In this paper we assumed that all low-density models (01, 02 and L) contain 
only compact lens objects with the same mass and radius, and compared their optical 
quantities. If we assume that particles in model 02 consist of compact lens objects 
and clouds with same mass but larger radius and that the number of compact lens 
objects is equal to that in models 01 and L, the values of optical quantities in 
model 02 will decrease as the radius of clouds increase, just as S(b) and S(c) in the 
Einstein-de Sitter model S, which were shown in the previous paper. 
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Fig. 3. The angular dependence of < >^^^. Solid and dotted lines denote behaviors for z = 1 
and 2, respectively. The separation angle 9 is in the unit of arcsec. 

For the statistical analysis we used 200 ray bundles reaching an observer in 
a single inhomogeneous model universe. This number of ray bundles may be too 
small to cover the influences from complicated inhomogeneities in all directions. For 
getting more robust statistical results it may be necessary to use more ray bundles 
and more model universes produced with random numbers. 

Here we touch recent observations of cosmological shear due to weak lensing 
and their relation to our results. Fort et al.tj) attempted the measurements of a 
coherent shear from foregroundjnass condensations in the fields of several luminous 
radio sources. Schneider et al.liil' determined the shear in the field (2 min x2 min) 
containing a radio source PKS 1508-05 with z = 1.2, and their result is that the 
shear is about 0.03 for the angular scale 1 min. If we interpret this value as an upper 
limit of the average shear in the general field due to weak lensing, we find from our 
result in Fig. 4 that Qq is ~ 0.4 at large. For its more confident estimate, more 
observational data of coherent shears in various fields are required. 
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Fig. 4. The angular dependence of < 7^ >^/^. Solid and dotted lines denote behavior for z = 1 
and 2, respectively. The separation angle 9 is in the unit of arcsec. 



Appendix A 

Basic Equations of Null Rays 



Corresponding to the line-element (2-10), we consider first the wave- vectors 
K'^ = {K^,K''') = dy^/dv, where v is an affine parameter. The null condition 
is 



C = -[Qq + {1-Qq- Ao)5 + Ao53]-ie6^(l + acP){Ky 



(A-l) 



where F is defined in Eg. ( 2-12 ). From the equation 

d dC _ dC 
dv'dIO' ~ dy^' 

we obtain 

di^o 3Qo + 2(1 - - Ao)e2(^-^o) ^ ^ 2G{T) 



(A-2) 



dv 



+ 



G{T) 



i\2 



1 



+ "^2 5r 



oy^ 
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Fig. 5. The z dependence of the angular diameter distance Da for a = 1. Models S(1,0), 01(0.2,0), 
02(0.4,0) and L(0.2,0.8) are denoted by solid, dotted, short dashed, long dashed lines, respec- 
tively. 



2cR-'e-'^cPG{T)J2{Kr/F'^=0, 



(A-3) 



H — 
2 dy' 



0, (A-4) 



where G{T) is defined by Eq. (|2-15| ). 

Transforming v to T, these equations lead to 



If 



+ 



3f?o + 2(1 -Qq- Ao)e^(^-^") 



+ 2 



(A-5) 



+ 4 - a—- if* - a- . 



+ {R^Ho/cf{l-no-Xo)F 



-1 
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1 



i 1^0 



0. 



(A-6) 



Moreover, if we use for 



K' = CR-^e-^K'/K^ 



(A-7) 



we obtain Eqs. (2-13) and (2-14) 
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Table I. Optical quantities in model Ol. 



z 


< K > 


< >i/^ 


< ^ > 


< >^'-' 


< 71 > 


< 7' 


1 


-0.0002 


0.0425 


1.0221 


0.1219 


-0.0028 


0.0426 


2 


-0.0054 


0.0796 


0.9961 


0.2210 


-0.0041 


0.0832 


3 


-0.0069 


0.1157 


1.0368 


0.6234 


0.0057 


0.1225 


4 


-0.0081 


0.1519 


1.0790 


0.5638 


0.0078 


0.1602 


5 


-0.0064 


0.1872 


1.1627 


1.1880 


0.0120 


0.1980 



Table II. Optical quantities in model 02. 



z 


< K > 


< ^ >^''' 


< M > 


< (1-m)' >^''' 


< 71 > 


< 7' >^''' 


1 


0.0008 


0.0471 


1.0126 


0.1305 


-0.0011 


0.0451 


2 


0.0071 


0.0842 


1.0517 


0.2755 


-0.0001 


0.0873 


3 


0.0111 


0.1103 


1.0980 


0.4956 


-0.0022 


0.1255 


4 


0.0136 


0.1380 


1.1602 


1.0742 


-0.0019 


0.1613 


5 


0.0147 


0.1680 


1.1926 


1.5755 


-0.0028 


0.1958 



Table III. Optical quantities in model L. 



z 


< K> 


< >^/'' 


< ^ > 


< >^'-' 


< 71 > 


<7' >V2 


1 


-0.0011 


0.0324 


1.0025 


0.0705 


0.0008 


0.0347 


2 


-0.0036 


0.0675 


1.0134 


0.1572 


0.0018 


0.0714 


3 


-0.0045 


0.0974 


1.0345 


0.2407 


0.0020 


0.1002 


4 


-0.0065 


0.1216 


1.0560 


0.3209 


0.0023 


0.1214 


5 


-0.0070 


0.1405 


1.0805 


0.4000 


0.0030 


0.1372 



